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I.  INTRODUCTION 


One  of  the  major  problems  facing  chemical  engineers 
involves  the  correlation,  extension,  and  prediction  of  phase 
equilibrium  data.     Except  for  systems  exhibiting  ideal  liquid 
phase  behavior,  this  problem  involves  an  understanding  of  the 
manner  in  which  activity  coefficients  vary  with  composition 
(and  temperatvire )  in  binary  and  multlcomponent  systems, 

Prausnitz  and  coworkers  [l]  have  presented  ample  evidence', 
to  show  that  activity  coefficients  in  systems  which  deviate 
only  moderately  from  ideal  behavior  can  be  well-represented  by 
the  Wilson  equation  [2],     These  workers  further  showed  that 
multicomponent  vapor-liquid  equilibrium  data  calculated  by 
means  of  the  Wilson  equation  with  constants  determined  from 
binary  data  were  in  agreement  with  experimental  measurements . 

Two  problems  arise  when  considering  partially  miscible 
systems s     (1)  Deviations  from  ideal  behavior  are  so  large  that 
it  is  difficult  to  fit  activity  coefficient  data  to  an  analyt- 
ical expression  of  moderate  complexity,    (2)  Sufficient  reliable 
equilibrium  data  for  partially  miscible  systems  are  not  avail- 
able so  that  the  behavior  of  these  systems  can  be  characterized. 

Renon  [3]  has  attacked  the  first  problem  by  deriving  a 
semi-theoretical  equation  based  on  Scott's  two-liquid  model 
and  a  nonrandom  mixture.     This  equation,  refered  to  as  the  NRTL 
(nonrandom,  two-liquid)  equation,  was  shown  to  be  capable  of 
representing  activity  coefficients  in  systems  exhibiting  large 
positive  deviations  from  ideal  behavior  when  tested  against 


the  few  available  experimental  data. 

The  problem        insufficient  experimental  data  for  partial- 
ly misclble  systc   -  has  received  little  attention  although 
these  systems  are        considerable  practical  interest.  Probably 
the  major  reason  for  the  neglect  of  these  systems  is  the  experi- 
mental difficulties  encountered. 

The  goals  of  this  study  weresj    (i)  to  develop  apparatus  and 
experimental  techniques  for  studying  these  difficult  systems 
and  (2)  to  obtain  experimental  data  so  that  the  behavior  of  par- 
tially misclble  systems  may  be  better  characterized. 

Systems  composed  of  Acetonitrile ,  n-Hexane  and  Cyclohexane 
were  chosen  for  study  because  there  is  a  relatively  large  degree 
of  mutual  solubility  and  because  the  relatively  high  volatility 
of  the  components  is  advantageous  for  the  equilibrium  measure- 
ments.    Isothermal  vapor-liquid  equilibrium  data  (including  the 
two-liquid  phase  region)  were  determined  for  the  n-hexane— 
acetonitrile  system  and  the  cyclohexane— acetonitrile  system. 
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II.  EXPERIMENTAL 

2.1  Apparatus 

The  equilibrium  total  vapor  pressure  measurements  were 
conducted  by  the  static  method.     The  schematic  arrangement  of 
the  apparatus  is  sho^vn  in  Figure  2-1,     The  apparatus  consisted 
of  a  thermostat  containing  an  equilibrium  cell  (1),  a  differen- 
tial pressure  transducer  (2)  and  a  thermometer  (3),  an  absolute 
pressure  gauge  (4),  a  differential  pressure  transducer  indi-  ' 
cator  (5),  a  pressure  regulating  bellows  (6)  and  a  vacuiim  system. 
Thermostat »     The  thermostat  was  an  insulated  l6"xl6"x3B"  box. 
On  the  front  side  of  the  box  there  was  a  10"x26"  acrylic  plastic 
pane  which  could  be  removed  for  servicing  the  equilibrium  cell 
during  shut-down.     Diiring  a  normal  run,  operation  of  the  valves 
and  stopcocks  located  inside  the  thermostat  were  done  by  hand 
through  the  isolating  gloves  hanging  on  the  side  walls  of  the 
box, 

A  6"  diameter  fan  circulated  the  air  upward  along  the  cen- 
ter and  downward  along  the  side  of  the  thermostat.     Around  the 
fan  were  two  nichrome  heaters,  one  of  them  was  regulated  manual- 
ly by  a  rheostat  and  the  other  was  automatically  controlled  by 
a  Ihermistep  Model  63  Temperature  Controller,     The  temperature 
variation  of  the  thermostat  was  less  than  ±  0.01°  C.  while 
maintained  at  ^5°  C.,  and  was  almost  at  a  uniform  temperature 
throughout  the  system. 

EquilibriuTTi  Cell;     The  equilibrium  cell  was  as  "shown  in  Figure 
2-2.     It  was  constructed  of  two  major  parts  —  the  vapor  sample 
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collector  and  the  equilibrium  chamber,  and  was  made  of  Pyrex 
glass.     The  vapor  sample  collector  consisted  of  two  sample  loops 
each  made  from  4  feet  of  size  C  capillary  tubing  and  sealed  to 
a  four-way,  twin-V  bore  stopcock.     The  loops  were  wound  together 
into    the  form  of  a  coil.     An  insulated  nichrome  wire  was  wound 
around  the  sample  loops  so  that  the  loops  could  be  heated  by 
passing  an  electric  current  through  the  wire. 

The  cell  was  connected  through  a  S  12/30  ground-glass-metal 
joint  to  the  differential  pressure  transducer  by  1/8"  O.D. 
copper  tubing.     In  the  equilibrium  chamber , there  were  two  liquid 
sampling  tubes  made  of  size  C  capillary  tubing.     These  liquid 
sampling  tubes  were  closed  by  Fischer  &  Porter  Lab-Crest  l^mm 
needle  valves  and  terminated  in  female  Luer  joints. 

The  vapor  sample  collector  and  the  equilibrium  chamber  were 
coupled  together  by  a  S  ^5/50  ground-glass  joint.     When  the 
determination  of  vapor  composition  was  not  required,  the  vapor 
sample  collector  was  replaced  by  a  cap.' 

The  equilibrium  cell  was  mounted  on  a  shaker  driven  by  a 
motor  located  on  the  top  of  the  box.     A  sketch  of  the  shaker 
mechanism  is  shown  in  Figure  2-3.     The  speed  of  the  motor  was 
regulated  by  a  Gerald  K.  Heller  Co. 's  Model  SXlOA  Motor  Contrdl- 
ler , 

Temperature  Measurement ;     A  Fisher  Cat.  No.  15-0^3  glass  thermo- 
meter with  0.1°  C.  subdivisions  was  used  to  measure  temperature. 
The  thermometer  was  crosschecked  with  several  other  similar 
thermometers  and  the  error  at  ^5°  C.  was  estimated  to  be  less 
than  +  0.05°  C. 


Fig.  2-3.  Shaker  of  Equilibrium  Cell 
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Pressure  Measurement:    A  pressure  transducer  and  an  Indicator 
manufactured  by  Pace  Engineering  Company  were  used  to  detect  the 
pressure  difference  between  the  equilibrium  cell  and  the  absolute 
pressure  gauge.     The  sensitivity  of  the  transducer  was  better 
than  ±  0.1  mm  Hg.     With  the  transducer  by-pass  valve,  V^,  closed, 
the  air  pressure  in  the  absolute  pressure  gauge  was  adjusted  by 
varying  the  volume  of  the  pressure  regulating  bellows  until  the 
reading  on  the  transducer  indicator  came  to  zero.     With  the 
transducer  indicating  zero  pressure  differential,  the  vapor  pres- 
sure in  the  equilibrium  cell  could  then  be  read  indirectly  from 
the  pressure  gauge.     The  absolute  pressure  gauge  was    a  Wallace 
&  Tiernan  Model  No.  FA129,  with  1  mm  Hg    subdivisions.     The  gauge 
had  been  calibrated  against  a  mercury  manometer-cathetometer 
system  and  required  no  correction. 

Vacuum  System;     The  vacum  pump  was  capable  of  pumping  down  to 
about  1  mm  Hg    absolute.     A  liquid  nitrogen  cooled  vapor  trap 
was  placed  ahead  the  pump, 
2,2     Operating  Procedure 

The  acetonitrile  used  in  this  work  was  Fisher  Scientific 
Company's  reagent  grade.     The  n-hexane,  cyclohexane,  and  benzene 
(used  as  a  solvent  in  the  calibration  of  the  gas  chromatograph) 
were  Phillips  Petroleum  Company's  "pure"  grade  reported  to  be 
99  mole  %  pure.     These  reagents  were  used  without  further  puri- 
fication because  chromatographic  analysis  revealed  no  signifi- 
cant impurities. 

Preparation  of  Equilibrium  Cell;     The  vapor  sample  collector 
was  filled  with  mercury  before  each  run.     The  device  used  for 
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filling  the  collector  with  mercury  is  illustrated  in  Figure  2-^. 

A  syringe  equipped    with  a  two-way  stopcock,  V;i_,  was  filled 
with  mercury  and  placed  on  the  terminal  "a"  of  the  vapor  sample 
collector.     With  the  stopcocks  turned  to  the  positions  as  illus- 
trated in  Figure  2->^,  the  terminal  "h^"  of  the  collector  v:as 
attached  to  a  vacuum  pump  and  one  half  of  the  apparatus,  in- 
cluding one  sample  loop,  was  evacuated.     Then  V2  was  closed 
opened,  and  the  mercury  was  introduced  into  the  loop  via  the 
syringe.     The  procedure  was  repeated  to  fill  the  other  half  of 
the  apparatus  with  mercury.     After  the  entire  apparatus  was 
filled  with  mercury,  the  syringe  was  removed  and  the  vacuum  pump 
was  disconnected.     With  the  four-way  stopcocks  remaining  in  the 
position  of  V3  in  Figure  2-4,  the  mercury  in  the  apparatus  was 
warmed  up  to  a  temperatiire  above  ^5°  C.  and  then  the  four-way 
stopcocks  were  rotated  to  seal  off  the  sample  loops  while  the 
mercury  was  warm. 

Before  charging  the  sepcimen    solution  into  the  equilibriim 
cell,  the  liquid  sampling  tubes  were  completely  purged  with  dry 
air  to  assure  that  no  solution  of  the  previous  rim  remained  in 
the  tubes.     If  the  specimen  solution  was  a  two-liquid  phase 
mixture,  about  50  ml  of  the  solution  would  be  charged,  and  if 
it  was  a  single  phase  mixture,  about  30  ml  of  the  solution  would 
be  charged.     The  vapor  sample  collector  and  the  equilibrim 
chamber  were  coupled  together  and  the  Joint  was  further  sealed 
by  applying  a  coat  of  Fisher  Sealit  over  the  seam  of  the  Joint. 
Derassing;;     After  the  apparatus  had  been  assembled,  the  foreign 
gases    present  in  the  vapor  space  as  well  as    dissolved  in  the 
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specimen  solution  were  removed  from  the  equilibrium  cell.  This 
vjas  accomplished  by  freezing  the  specimen  solution  with  liquid 
nitrogen  until  it  had  become  solid,  then  evacuating  the  cell 
through  the  vacuum  system.     During  the  evacuation,  the  transducer 
by-pass  valve,  V]_,  and  valve  V2  in  Figure  2-1  were  open.  After 
the  evacuation  of  the  cell  was  complete,  valve  V2  was  closed  and  ■ 
the  frozen  specimen  solution  allowed  to  melt.     When  this  process  ' 
was  repeated  three  times,  the  equilibrium  cell  together  with  the 
specimen  solution  was  effectively  degassed,  and  no  bubble  forma- 
tion during  normal  freezing    was  observed.     Before  evacuating 
the  system,  care  was  exercised  to  assure  that  the  two-way  stop- 
cocks on  the  vapor  sample  collector  were  closed.     The  melting 
of  the  frozen  specimen    was  accelerated  with  the  aid  of  a  port- 
able hot  air  blower.     At  the  end  of  the  degassing  process  the 
transducer  by-pass  valve,  Yi,  was  closed  instead  of  valve  V2, 
and  the  four-way  stopcocks  of  the  vapor  sample  collector  were 
turned  to  seal  off  both  the  sample  loops  and  the  mercury  bulbs. 
Determination  of  Equlllbrlim  and  Collections  of  Samples;  After 
the  temperature  of  the  thermostat  had  been  brought  to  ^5°  C, 
the  equlllbriiom  cell  was  then  shaken  frequently  to  agitate  its 
contents.     In  the  mean  time  the  air  pressure  In  the  absolute 
pressure  gauge  was  adjusted  so  as  to  maintain  the  pressure 
difference  between  the  equilibrium  cell  and  the  absolute  pressure 
gauge  at  a  minimum  until  no  further  change  in  vapor  pressure  in 
the  equllibrivim  cell  was  observed,  then  the  transducer  indicator 
was  adjusted  to  zero  and  the  equilibrium  total  vapor  pressure 
of  the  mixture  was-  determined. 
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Before  introducing  vapor  samples  into  the  sample  loops,  the 
sample  inlet  tubes  were  purged  by  dropping  out  about  2/3  of  the 
mercTJiry  stored  in  the  mercury  bulbs  with  the  four-way  stopcocks 
turned  to  the  position  of  V]_  in  Figure  2-8 (A)  and  by  opening  the 
two-way  stopcocks.     After  the  sample  inlet  tubes  had  been  purged, 
the  four-way  stopcocks  were  turned  to  the  position  of  V]_  in 
Figure  2-8(B),  then  the  vapor  sample  was  introduced  into  the 
sample  loops  while  the  mercury  in  the  loops  dropped  out.  After 
the  vapor  sample  had  been  collected,  the  four-ivay  stopcocks 
were  then  rotated  to  seal  off  the  sample  loops. 

For  collecting  liquid  samples,  the  equilibrium  cell  vras 
first  pressurized  to  about  10  to  20  mm  Hg    above  the  atmospheric 
pressure  by  introducing  dry  air  from  outside,  then  the  liquid 
samples  were  pushed  out  of  the  cell  through  the  liquid  sampling 
tubes  and  introduced  into  l/l6"  CD.  indium  tubing  attached  to 
the  sampling  terminals  as  shown  in  Figure  2-5.     The  liquid 
samples  were  then  sealed  into  the  indium  tube  capsules  of  about 

5  mm  long,  using  a  special  crimping  tool  manufactured  by  Sargent 

6  Company.     Usioally  10  to  15  capsules  for  each  liquid  sample 
v^ere  prepared  for  analysis. 

During  the  collection  of  samples,  the  temperature  of  the 
thermostat  was  maintained  at  ^5°  C.     Before  the  vapor  sample 
collector  was  removed  from  the  thermostat  for  analysis,  the 
sample  loops  were  warmed  up  to  above  ^5°  C.  and  were  kept  warm 
during  the  analysis  to  prevent  vapor  condensation. 
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Needle  Vaiva 


Fig.  2-5.  Collection  of  liquid 
samples  . 
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2.3  Analysis 

The  vapor  and  liquid  samples  were  analyzed  chromatographi- 
cally  to  determine  their  compositions.     A  Varian  Aerograph 
Model  600  D  Gas  Chromatograph,  associated  v;ith  a  Sargent  Model 
SR  Recorder  and  a  Disc  Integrator  Model  204,  was  used  for  the 
analysis. 

The  chromatographic  column  was  a  5  ft.xl/8  in.  O.D.  stain- 
less steel  tube  packed  with  20^  DECS  on  6O/8O  Chrom  W  supplied 
by  the  Varian  Aerograph  Company.     The  detector  was  a  flame 
ionization  detector,  and  nitrogen  was  the  carrier  gas.  The 
temperature  of  the  column  was  maintained  at  60°  C.  and  the  liquid 
sample  size  was  approximately  1  microliter. 

The  ratios  of  the  chromatogram  peak  areas  for  n-hexane  and 
acetonitrile  as  well  as  cyclohexane  and  acetonitrile  were  found 
to  be  sensitive  to  changes  in  liquid  sample  size,  especially 
when  the  sample  sizes  were  smaller  than  0,8  microliter,  however, 
no  significant  effect  was  observed  if  sample  sizes  were  within 
the  range  of  0.8  to  I.5  microliters. 

The  correction  factors  which  converted  the  area  ratios  of 
the  chromatographic  peaks  to  the  weight  ratios  of  the  components 
for  n-hexane — acetonitrile  and  cyclohexane — acetonitrile  systems 
viere  found  to  vary  with  the  composition  of  the  mixture.  Cali- 
bration curves  for  both  systems  were-  constructed  by  use  of  liquid 
mixtures  of  known  composition  prepared  by  weighing.  To  the  partial 
ly  miscible  mixtures,  benzene  was  added  as  a  solvent  so  as  to 
convert  the  mixtures  into  single  phases.     The  calibration  curves 
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are  shown  in  Figures  2-6  and.  2-?,  The  standard  relative  deviation 
°rel  •  chromatographic  peak  area  ratios  was  normally  atiout 

3  to  3.5^.     Therefore,  an  arithematlc  average  was  taken  of  8  to 
10  measurements  to  determine  a  correction  factor  or  a  composition. 

To  introduce  the  vapor  sample  into  the  chromatograph ,  the 
terminals  "a"  and  "b-j_"  of  the  vapor  sample  collector  as  shown  in 
Figure  2-8  were  tied  into  the  carrier  gas  stream  with  Teflon 
tubing.     With  the  stopcocks  in  the  positions  shown  in  Figure 
2-8(A),  carrier  gas  was  swept  through  b2^C]_a  to  purge  the  sample 
lines.     The  four-way  stopcock  Y-^  then  momentarily  turned  to 

the  position  shown  in  Figure  2-8(B),  and  a  plug  of  gas  in  the 
sample  loop  was  swept  onto  the  chromatographic  column.     The  amount 
of  the  sample  gas  swept  onto  the  column  depended ■ upon  the  length 
of  time  stopcock  V^^  remained  in  the  position  shown  in  Figure 
2-8(B).     It  was  usually  possible  to  obtain  four  or  five  measure- 
ments from  one  sample  loop. 

To  introduce  the  liquid  sample  into  the  chromatograph,  the 
indium  capsule  containing  the  liquid  sample  was  melted  in  a 
furnace  swept  by  carrier  gas  and  the  sample  was  released  and  ■ 
vaporized.     The  apparatus  used  for  introducing  the  liquid  sample 
into  the  chromatograph  was  an  Aerograph  Inductor  Model  No.  695 
manufactured  by  Wilkens  Instrument  and  Reasearch  Inc.  The 
temperature  of  the  furnace  was  maintained  at  350°  C. 
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Fig.  2-6.  Gas  Chromatograph  Calibration  Curve 
System:  n-Hexane — Acetonitrile 
Columns  5'xl/8"  SS  20%  DECS  on  60/80 

Chrora  W 
Column  Temperature:  60°  C. 
Sample  Size:  1  microliter  (liquid). 
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Fig.  2-7,  Gas  Chromatograph  Calibration  Curve 
System:  Cyclohexane — Acetonltrile 
Column:  5'xl/8"  SS  20^  DECS  on  60/80 

Chrom  W 
Column  Temperature:  60°  C. 
Sample  Size:  1  microliter  (liquid). 
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2.4  Data 

The  experimentally  determined  equilibrium  total  vapor 
pressure-liquid  composition  data  are  listed  in  Tables  2-1  and 
2-2,  and  plotted  in  Figures  3-5(a),(b)  and  3-6(a),(b). 

For  binary  systems,  vapor  phase  compositions  can  be  calcu- 
lated by  rigorous  thermodynamic  equations  from  the  total  pres- 
sure-liquid composition  data  measured  at  a  constant  temperature 
and  hence  direct  determination  of  vapor  compositions  is  un- 
necessary.    The  apparatus  and  experimental  techniques  developed 
f ordetermining  mult i component  vapor  compositions  at  low  pressure 
were  tested  for  their  acciiracy  and  precision  on  the  n-hexane  — 
acetonitrile  system.     These  vapor  composition  data  are  also 
listed  in  Table  2-1,  and  will  be  compared  with  those  calculated 
from  the  thermodynamic  equations  in  Chapter  3. 
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Table  2-1.  Total  Vapor  Pressure-Liquid  Composition 
-Vapor  Composition  Measurements 

System  I     n-Hexane — Acetonitrlle  at  ^5°  C, 

P  XI  yi 

mm  Hg.abs,  Mole  Fraction  Mole  Fraction 

n-Hexane  in  Liquid      n-Hexane  In  Vapor 

213.0  e.OOOO  .000 

291.2  0.0135  .312 

300.3  •  '  0.0159  ■  .386 
31^.2  0.0203        .  .^19 

357.0  0.032^  •  .^55 
381. if  0.0397  .^87 
403.2  0.0^^5 

407.2  0.0502  -  .492 

420.1  0.0550        .        .  .507 

443.2  0.0638  .567 
459.0  0.0686  .■  .584 

462.0  ■  .  -0.0710 

482.1  0.0819 

493.0  •        0.0910  .598 
500.8  •  0.0945 

520.5  0.121       _____  . 

j;ii:=='  .643 

520.5  0.8525   

5I8.7  0.8780  ■ 

509.1  0.9154  .     ■       ■•  ■ 
497.0  .  0.9305 

491.6  0.9416  •  V 

473.3  0.9539 
446.0  •  0.9706 

343.3  1.00.00 
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Table  2-2  .  Total  Vapor  Pressure-Liquid  Composition 
Measurements 

System:     Cyclohexane— Acetonitrile  at  ^5°  C. 


P  XI 
mm  Hg.  Absolute  Mole  Fraction 

Cyclohexane  in  Liquid 


213.0 

0.0000 

n  o  O  A 

239.0 

0 . Olio 

0  CO.  1 

A     A  T  A  A 

0 . 0190 

A     A  0  J 1  ^ 

0.0290 

0.0339 

309. 7 

0 .  O^lll 

TOT  h 

A    A    A  n 

0  •  0501 

33j»^ 

A     A  CT  ^ 

A     A/^  A  A 

365.0 

0.0715 

37^.8 

0.0755 

370.3 

0.0760 

390.0 

O.O867 

396.6 

0.0990 

^06.3 

0.106 

ii-08.5 

0.112 

^09.0 

0.120 

409.0 

0.8877 

403.0 

0.9081 

395.0 

0.9195 

388.8  . 

■    •  0.9271 

377.3  . 

0.9418 

344.6  ■ 

0.9644 

328.0 

0.9686 

325.2 

0.9793 

233.7 

1.0000 
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III.     DETERMINATION  OF  COMPLETE  VAPOR-LIQUID  EQUILIBRIUM 
DATA  OF  THE  BINARY  SYSTEMS  FROM  TOTAL  PRESSURE- 
LIQUID  COMPOSITION  MEASUREMENTS 

3 . 1    Theoretical  Backgrround 

Determination  of  vapor-liquid  equilibrium  data  from  total 
pressure-liquid  composition  measurements  can  be  achieved  either 
by  (1)  direct  step-wise  integration  of  the  Gibbs-Duhem  equation 
as  illustrated  by  Othraer  et.  al.  [^],   or  (2)  evaluation  of  the 
constants  of  equations  representing  the  liquid  activity  coef- 
ficients by  the  method  of  least  squares  as  suggested  by  Barker 
[5]. 

Direct  Evaluation  of  Equilibrium  Vapor  Composition;  From  thermo- 
dynamic reasoning,  the  general  form  of  the  Glbbs-Duhem  equation 

SdT  -  VdP  +  £n3_dHj_  =  0  (3-1) 
implies  a  relation  between  simultaneous  changes  of  temperature,  ■ 
pressure  and  the  chemical  potentials  of  the  individi;ial  consti- 
tuents of  a  single  phase  mixture.     At  constant  temperature, 
equation  (3-1)  reduces  to 

Sn^djii  =  VdP  (3-2). 
This  eqmtion  may,  of  course,  be  applied  to  any  phase,  however, 
here  it  is  applied  specifically  to  the  liquid  phase.  Therefore, 
V  is  the  volume  of  the  liquid  solution. 

For  a  system  x^rith  liquid  and  vapor  phases  in  equilibriiim, 

the  criterion  of  equilibrium  is 

liquid  vapor 
;^i  =  ?i  (3-3). 

Further,  if  the  vapor  phase  is  a  perfect  gas  mixture,  the  chemical 
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potentials  may  be  expressed  as 

=  }x°  +  RT  lnpj_  (3-^), 
where        is  a  function  only  of  temperature.    At  constant  tempera- 
ture , 

djx^  =  RT  dlnpi  ,         /  •  ■  (3-5). 

Substituting  (3-5)  in  (3-2)  gives 

RTIn^dlnp^  =  VdP  (3-6).  ' 

Equation  (3-6)  gives  a  relationship  between  changes  in  the  total 
pressure  P  and  the  partial  pressures  Pj_. 

Dividing  equation  (3-6)  through  by  the  total  niimber  of  moles, 
In-j_,  in  the  condensed  phase,  we  obtain 

Xx^dlnpi  =  VidP/RT  ,  (3-7). 

Let  RT  =  PVg,  then  equation  (3-7)  may  be  rewritten  as 

Ixidlnpi  =  (vi/Vg)dlnP  (3-8). 

For  a  binary  system  with  components  1  and  2,  the  partial 
pressures  are  p^  =  y^P  and  P2  =  (1  -  yi)P,  respectively. 
Substituting  these  relations  in  (3-8),  gives  \ 

Xidln(y]^P)  +  (l-xi)dln[(l-y-L)p]  =  (v;i_/Vg)dlnP  (3-9). 
And,  on  rearranging  this  equation  we  obtain 

(dlnP)T  =   li-  11  _  dyi  (3-10). 

yi(i  -  yi)Ci  -  (v/Vg)] 

Equation  (3-10)  gives  the  relationship  between  the  total 
vapor  pressure  measurements  and  the  compositions  of  liquid  and 
vapor  phases  of  a  binary  system  in  equilibrlim  at  constant  tem- 
perature.    At  low  pressure  the  term  [l  -  (V]^/Vg)3  may  be  consid- 
ered unity  with  little  error.     Then,  the  equation  simplifies  to 


2k 

(dlnP)T,  =         "  dy,  (3-11). 

^      Yld  -  yi)  ^ 

The  application  of  equation  (3-11)  makes  possible  the  rig- 
orous thermodynamic  calculation  of  vapor  composition  data  from 
total  vapor  pressure-liquid  composition  jiata. 

Evaluation  of  Constants  of  Equations  Representing;  Activity  Coef- 
ficients ;     For  a  binary  liquid  system  in  equillbriim  with  an 
ideal  vapor  phase,  the  relationship  between  liquid  composition 
and  total  vapor  pressure  is 

P  =  XiYipf  +  (1  -  xi)Y2P2°  (3-12). 
It  is  assumed  that  the  activity  coefficients  y^  and  Y2  equation 
(3-12)  may  be  represented  by  a  particular  activity  coefficient 
equation.     The  problem  is  therefore  to  determine  the  constants 
in  the  activity  coefficient  equation  by  a  process  of  successive 
approximations  so  that  the  total  vapor  pressure  values  calculated 
from  equation  (3-12)  agree  with  the  experimentally  observed  values. 

The  procedure  for  earring  out  these   calculations  has  been 
detailed  by  Barker  [53*  and  a  complete  computer  program    has  been 
published  by  Prausnitz  et.  al.  [l]. 
3.2    Calculation  and  Results 

Ster?-wise  InteG:ration  of  Gibbs-Duhem  Equation?     In  the  calculation 
of  equilibrium  vapor  composition  data  from  total  vapor  pressure- 
liquid  composition  data  for  the  n-hexane — acetonitrlle  system  and 
the  cyclohexane — acetonitrlle  system,  the  step-wise  integration 
procedure  was  carried  out  by  using  eq-uatlon  (3-11)  in  the  form  of 
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61nP  =  lllJlillL  6y-,  (3-13). 

The  estimation  of  the  vapor  phase  fugacity  coefficients 
of  both  n-hexane — acetonitrile  and  cyclohexane — acetonitrile 
systems,  according  to  the  method  suggested  by  Prausnitz  et.  al. 
[l]  and  using  the  second  virial  coefficient  of  acetonitrile 
reported  by  Prausnitz  and  Carter  [6],   justified  the  ideal 
vapor  phase  assiamption  made  in  the  derivation  of  equation  (3-13)  t 

An  increment  of  6xi  =  0.005  mole  fraction  was  employed  in 
the  step-wise  calculation.     Values  of  P  and  x  needed  for  the 
calculation  were  taken  from  smooth  curves  drawn  through  the 
experimental  data  points.     The  integration  was  started  at 
x^^*^^  =  0  where  y^^^^^  =  0  and  P^°^  =  P2°.     For  the  first  increment  we 
have 

61nP  =  InP^  -  lnP^°^ 

A  CD 
6yi  =  yi"^ 

and , 

CD  _  61nP  ^  xi^^ 
1  +  61nP 

The  equilibrium  vapor  composition  at  the  end  of  the  first  incre- 
ment, y-^CD,  was  thus  obtained.     At  the  end  of  the  second  increment, 
the  equilibrium  vapor  composition  was  obtained  by  trial  and  error 
calculation  for  y-^^^    from  equation  (3-13)  with  x^^^)  _  o.Ol  and 

61nP  =  InP^^  -  InP^ 

6y^  =  y^^)  -  y ,CD. 

Subsequent  integrations  were  carried  out  in  the  same  manner  until 
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Xq_  reached  the  solubility  limits. 

Equilibrium  vapor  compositions  on  the  other  side  of  the 
miscibility  gap  were  obtained  by  starting  the  integration  from 
X]_  =  1.0.     The  results  are  listed  in  Tables  3-1  and  3-2,  and 
plotted  in  Figiires  3-1  and  3-2.     The  equilibrium  vapor  composi- 
tion should  remain  constant  over  the  two-liquid  phase  region 
■hence  this  composition  obtained  from  the  calculations  begining 
at  x-]_  =  0  should  match  that  obtained  from  the  calculations  begin- 
ning at  X]_  =  1,0.     Reasonable  matchings  were  obtained  for  both 
systems.     The  mismatch  in  the  n-hexane — acetonitrile  system  was 
only  0.016  mole  fraction,  while  in  the  cyclohexane — acetonitrile 
system  it  was  0.02?  mole  fraction. 

The  activity  coefficients  of  the  components  in  the  liquid 
phase  at  ^5°  C.  were  calculated  from  the  equilibrium  data  using 
following  relationships : 

Yi  =  Pyi/Pi°xi 

(3-15). 

Y2  =  Py2/P2°X2 

The  values  obtained  are  also  listed  in  Tables  3-1  and  3-2,  and 
plotted  in  Figures  3-3  and  3-4. 

Evaluation  of  Constants  in  Activity  Coefficient  Equations:  The 
constants  in  the  following  activity  coefficient  eqixations  were 
evaltiated  using  the  IBM  36O  computer    with  a  least  squares  program 
similar  to  that  of  Prausnitz  et.  al.  [l]. 
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Table  3-1.  Calculated  Equilibrium  Vapor  Compositions 
and  Activity  Coefficients  by  Step-v;ise 
Integration  of  Gibbs-Duhem  Eqviation 

System:     n-Hexane ( 1 ) — Acetonitrile (2 ) 
at  45°  C. 


1 

P 

1 

'  1 

'  2 

0.000 

213.0 

0.000 

0.005 

2^7 

0.134 

0.010 

273 

0.216 

17.2 

1.02 

0.015 

295 

0.275 

0.020 

315 

0.322  " 

14.8 

1.02 

0.025 

332 

0.358 

0.030 

3^9 

0.391 

.  13.3 

1.03 

0.035 

■  365 

0.419 

0.0^0 

380 

0 . 444 

12.3 

1.03 

0.0^5 

394 

0.465 

0.050 

409 

0.486 

11.6 

1.04 

u .  0^0 

0.060 

463 

0.521 

11.0 

1.04 

0  06'^ 

\J  %  \J  ^  _J 

44  R 

0.070 

459 

0.548 

10.5 

1.05 

0.075 

470 

0.560 

0.080 

479 

0. 570 

9.94 

1.06 

0.085 

488 

0.579 

0.090 

495 

0.586 

9.39 

1.06 

0.095 

502 

0.593 

0.100 

508 

0.599 

8.86 

1.06 

0.105 

513 

0 .  6o4 

0.110 

516 

**■  ^ 

0 . 607 

8.29 

1.07 

0.115 

519 

0 . 609 

0.121 

520,  5 

0.610 

7  64 

1.08 

id 

Xt 

1 

■^2 

•^1 

0.000 

1.000 

343.3 

0.000 

1.000 

0.005 

0.995 

365 

0.062 

0.938 

0.010 

0.990 

386 

0.116 

0.884 

1.00 

21,0 

0.015 

0.985 

4o4 

0.158 

0.842 

0.020 

0.980 

421 

0.195 

O.8O5 

1.01 

19.3 

0.025 

0.975 

435- 

0.224 

0.776 

0.030 

0.970 

44? 

0.247 

0.753 

1.01 

17.3 

0.035 

0.965 

457 

0.266 

0.734 

0 . 040 

0.960 

466 

0.282 

0.718 

1.02 

15.4 

0.0^5 

0.955 

474 

0.296 

0.704 

0.050 

0.950 

481 

O.3O8 

0.692 

1.02 

13.9 

0.055 

0.9^5 

487 

O.3I8 

0.682 

0.060 

0.9^0 

492 

0.326 

0.674 

1.03 

12.6 

0.070 

0.930 

500 

0 . 340 

0.660 

1.03 

11.4 

0.080 

0.920 

507 

0.351 

0.649 

1.04 

10.5 

0.090 

0.910 

512 

0.359 

0.641 

1.05 

9.56 
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rable  3-1.  Cont'd. 


^2 

^1 

P 

yi 

Yl 

Y2 

0.100 

0.900 

514 

0.362 

0.638 

1.06 

8.74 

0.110 

0.890 

517 

0.367 

0.633 

1.07 

8.10 

0.120 

0.880 

519 

0.371 

0.629 

1.08 

7.53 

0.130 

0.870 

520 

0.373 

0.627 

1.08 

7.01 

0.1475 

0.8525 

520.5 

0.374 

0.626 

1.11 

6.20 

Table  3-2.  Calculated  Equilibrium  Vapor  Compositions 
and  Activity  Coefficients  by  Step-wise 
Integration  of  Gibbs-Duhem  Equation 

System »     Cyclohexane (1 ) — Acetonitrlle (2 ) 
at  45°  C. 


"^1 

P 

yi 

Yi 

Y2 

0.000 

213.0 

0.000 

0.005 

225 

0.057 

0.010 

237 

0.108 

11.0 

1.00 

0.015 

249 

0.154 

0.020 

261 

0.196 

11.0 

1.01 

0.025 

273 

0.235 

0.030 

284 

0.268 

10,9 

1.01 

0.035 

295 

0.298 

0.040 

306 

0.326 

10.7 

1.01 

0.045 

317 

0.352 

0.050 

328 

0.376 

10.6 

1.01 

0.055 

338 

0.397 

0.060 

348 

0.417 

10.4 

1.01 

0.065 

357 

0.434 

0.070 

365 

0.448 

10.0 

1.02 

0.075 

373 

0.462 

0.080 

380 

0.474 

9.6 

1.02 

0.085 

387 

0.485 

0.090 

393 

0.494 

9.2 

1.03 

0.100 

402 

O.508 

8.7 

1.03 

0.110 

407 

0.516 

8.2 

1.04 

0.120 

409 

0.519 

7.6 

1.05 

Table  3-2.  Cont'd. 


c 

^1 

P 

c 

X 

X 

c 

0.000 

1.000 

233.7 

0.000 

1.000 

0.005 

0.995 

256 

0.088 

0.912 

0.010 

0.990 

278 

0.162 

0.838 

1.01 

21.2 

0.015 

O.9B5 

296 

0.215 

0.785 

0.020 

0.980 

310 

0.253 

0.747 

1.01 

18.4 

0.025 

0.975 

322 

0.283 

0.717 

16.3 

0.030 

0.970 

33^ 

0.311 

0.689 

1.02 

0.035 

0.965 

34^ 

0.333 

0.657 

0.0^0 

0.960 

352 

0.350 

0.650 

1.02 

14.5 

0.0^5 

0.955 

360 

0.366 

0.634 

0.050 

0.950 

367 

O.38O 

0.620 

1.03 

13.1 

0.055 

0.9^5 

37^ 

0.393 

0.607 

0.060 

0.9^0 

379 

0.^^02 

0.598 

1.03 

11.9 

0.065 

0.935 

38^ 

0.411 

0.589 

0.070 

0.930 

388 

0.418 

0.582 

1.04 

10.9 

0.075 

0.925 

393 

0.427 

0.573 

0.080 

0.920 

396 

0.432 

0.568 

1.05 

10.0 

0.085 

0.915 

399 

0.437 

0.563 

0.090 

0.910 

^02 

0.442 

0.558 

1.05 

9.27 

0.100 

0.900 

i^06 

0.449 

0.551 

1.06 

8.86 

0.113 

0.887 

409.0 

0.454 

0.546 

1.08 

7.72 
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(1 )  Three-parameter  T'Tar.yules  Equation 
InYi  =        [a  +  2x3_(B  -  A  -  D)  +  3X]_^  D] 
lnY2  =  +  2X2(A  -  B  -  D)  +  3x3^  D] 

( 2 )  Tvfo-parameter  Marp:ules  Equation 
InYi  =        [a  +  2xi(B  -  a)] 
lnY2  =        Lb  +  2x2 (A  -  B)] 

( 3 )  van  Laar  Equation 


-.2 

InYi  =  +  (A/B)xi/X2] 


/ 


lnY2  =  3/[l  +  (B/A)x2/xir 

(4)    NRTL  (nonrandom.  two-llquld)  Equation 

InYn  -  Xo2  /t  exp(-2Ai2T2i) 

inYi  =  ^2  \  ^217  — — - — T^2 

I       [xi  +  X2exp(-Ai2T2l)J 

exp(-Ai2Ti2) 


+  T12 


[xg  +  Xiexp(-Ai2Ti2)] 


2( 


lnY2  =  xi^ 


,  exp(-2Ai2Ti2) 
.  ]^2  '  2 

[x2  +  xiexp(-Ai2Ti2)] 

+  T21-  exp(-Ai2T2i) 


(3-16) 


(3-17) 


(3-18) 


(3-19) 


[xi  +  X2exp(-Ai2T2l)]  I 

The  constants  determined  are  listed  in  Table  3-3. 

The  equilibrium  total  vapor  pressures,  equilibrium  vapor 
compositions  and  liquid  phase  activity  coefficients  obtained  from 
the  fit  of  the  two-parameter  and  three-parameter  Kargules  equa- 
tions are  listed  in  Tables  3-4  and  3-5,  and  plotted  in  Figures 
3-5(a),(b),  3-6(a),(b),  3-1,  3-2,  3-3  and  3-^,  respectively,  for 
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comparison  v;ith  the  experimental  data  and  those  calculated  from 
the  step-v;ise  integration  of  the  Gibbs-Diihem  eqviation.  The 
results  obtained  from  the  fit  of  the  van  Laar  equation  were 
essentially  the  same  as  those  obtained  from  the  fit  of  the  two- 
parameter  Marsules  equation,  and  those  obtained  from  the  fit  of 
the  NRTL  equation  were  similar  to  those  obtained  from  the  fit  of 
the  three-parameter  Margules  equation  if  the  constant  k\2  T'^as 
not  fixed;  but  were  similar  to  those  obtained  from  the  fit  of 
the  two-parameter  Margules  equation  ^f  the  constant  A12  was 
fixed  at  0.200. 


Table  3-3,  Constants  of  activity  coefficient 
equations  evalxxated. 


n-Hexane — 

Acetonitrile 


Cyclohexane— 
Acetonitriles 


2-Parameter 

y.argules 
Equation 


A  =  1.1853 
B  =  1.2152 


A 
B 


1.1299 
1.2271 


3-Parameter 
Kar gules 
Eqiiation 


A  =  1,2655 
B  =  1.3217 
D  =  ,68263 


A  =  1,1587 
B  =  1.2570 
D  =  ,22239 


van  Laar 
Sqiiation 


A  =  I.I85I 
B  =  1.2158 


A  =  1.1307 
B  =1.2286 


NRTL 
Equation 


Ai2=  .37201 
Ti2=  2.1336 
T2i=  1.9528 


Ai2=  .26588 
Ti2=  1.8976 
T21=  1.52^5 


NRTL 
Equation 


Ai2=  0.2000 
Ti2=  1.7170 
T2i=  1.5^59 


A]_2=  0.2000 
Ti2=  1.8305 
T2i=  1.3672 
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3-l»  Equilibrium  composition  values  for  the  system 
n-Kexane-Acetonltrll©  at  ^50  C. 
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3-2.  Equilibrliim  compos  it  Ion  values  for  the  system 
Cyclohexane-Acetonitrlle  at  45^  C, 
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3.  Activity  Coefficients,  System  n-Hexane(l) — 
Acetoni'trile(2)  at  45°  c. 
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  Gibbs-Duhera  Equation 

 2-paranieter  Margules 

Equation 

 3-parameter  Margules 

Equation 
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x^,  mole  fraction  cyclohexane  in  liquid 

Fig.  3-^.  Activity  Coefficients.  System  Cyclohexane ( 1 )- 
Acetonitrile(2)  at  45°  C. 
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X,  mole  fraction  n-hexane  in  liquid 
Pift.   3-5£i«  Total  vapor  pressure  vs.  liquid  composition, 
system  n-Hexane-Acetonitrlle  at  ^5°  C. 
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X,  mole  fraction  n-hexane  in  liquid 
3-5'i3.  Total  vapor  pressure  vs.  liquid,  composition, 
system  n-Hexane-Acetonitrile    at  45©  C, 
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X,  mole  fraction  cyclohexane  in  liquid 
is:.   3-6a.     Total  vapor  pressure  vs.  liquid  composition, 
system  Cyclohexane-Acetonitrle  at  ^5°  c. 
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X,  mole  fraction  cyclohexane  in  liquid 
Fig.  3-ob.  Total  vapor  pressure  vs.  liquid  composition, 
system  Cyclohexane-Acetonitrile  at  ^50  c. 
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IV.     DISCUSSION  AND  CONCLUSIONS 

The  experimental  total  vapor  pressure-liquid  composition 
data  obtained  in  this  work  seemed  to  be  of  good  quality  because 
smooth  curves  were  obtained  on  their  plots  and  good  matchings  in 
the  vapor  compositions  over  the  two-liquid  phases  were  obtained 
with  the  direct  Gibbs-Duhem  equation  calculations.     The  fact 
that  the  activity  coefficient  equations  do  not  fit  the  total 
pressure-liquid  composition  data  as  well  as  expected  could  mean 
only  that  these  equations  are  simply  not  capable  of  representing 
activity  coefficients  sufficiently  well  in  systems  with  large 
positive  deviations  such  as  these. 

The  experimentally  determined  vapor  compositions  of  the 
n-hexane — acetonitrile  system  were  only  in  fair  agreement  with 
those  calculated  from  the  total  pressure-liquid  composition 
measurements  by  the  step-wise  integration  of  the  Gibbs-Duhem 
equation.     This  is  believed  to  be  due  to  errors  in  analysing  the 
vapor  samples.     The  fact  that  the  area  ratios  of  the  chromato- 
graphic peaks  were  sensitive  to  the  sample  size  and  the  correction 
factors  which  converted  the  area  ratios  to  the  weight  ratios  of 
the  components  were  sensitive  to  the  composition  of  the  mixture, 
for  both  n-hexane — acetonitrile  and  cyclohexane — acetonitrile 
systems,  indicated  the  possibility  of  partial  adsorption  of  the 
components  in  the  gas  chromatograph.     The  technique  used  in  intro- 
ducing vapor  samples  into  the  gas  chromatograph  could  not  ade- 
quately reproduce  the  corresponding  sample  sizes  used  in  the 
construction  of  the  calibration  curves.     Because  the  extent  of  the 


^3 


adsorption  effect  could  be  dependent  on  sample  size,  the  inability 
to  reproduce  the  sample  size  could  cause  errors  in  the  analyses. 
Improvement  of  the  technique  for  vapor  sample  introduction 
might  possibly  improve  the  accuracy  of  the  vapor  composition 
determination;  however,  the  selection  of  a  chromatographic  colximn 
which  will  not  adsorb  the  components  of  the  sample  introduced 
would  not  only  improve  the  accuracy  of  the  vapor  composition 
determinations,  but  also  improve  the  precision  of  the  chromato- 
graphic measurements. 
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NOMENCLATURE 


A  Constant  in  activity  coefficient  equation 

B  Constant  in  activity  coefficient  equation 

D  Constant  in  activity  coefficient  eqmtion 

n  Number  of  moles  of  component;  Number  of  data  points 

P  Total  vapor  pressure  of  system 

p  Partial  vapor  pressure  of  component 

R  Gas  constant 

5  Entropy  of  system 

T  Absolute  temperature;  Constant  in  NRTL  equation 

V  Total  volume 

Vg  Voliame  per  mole  of  vapor  phase 

v^  Volume  per  mole  of  condensed  phase 

X  Mole  fraction  in  condensed  phase 

y  Mole  fraction  in  vapor  phase 

Y  Liquid  phase  activity  coefficient 

6  Sign  indicating  excess  of  final  over  initial  value 
p.  Chemical  potential  . 

^rel.  Standard  relative  deviation 

Z  Summation  operator  sign 
Subscripts 

♦ 

1  Component  1 

2  Component  2 

i  1^^  component 
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Apparatus  and  experimental  techniques  were  developed  for 
studying  phase  equilibrium  of  partially  miscible  systems.  Vapor- 
liquid  equilibrium  data  (including  the  two-liquid  phase  region) 
were  obtained  for  the  binary  systems  n-hexane — acetonitrile  and 
cyclohexane — acetonitrile  at  ^5°  C. 

The  equilibrium  total  vapor  pressure  measurements  vxere 
conducted  by  the  static  method  using  a  nevjly  designed  equi- 
libriiim  cell,  and  the  corresponding  equilibrium  vapor  and  liquid 
compositions  were  determined  by         gas  chromatographic  analysis. 
The  equilibrium  vapor  compositions  and  the  liquid  phase  activity 
coefficients  were  also  calculated  from  the  total  vapor  pressin'e- 
liquid  composition  measurements  by  (1)  the  step-wise  integration 
of  the  Gibbs-Duhem  equation  and  (2)  the  non-linear  least  squares 
fit  of  the  activity  coefficient  equations  including  the  three- 
parameter  and  two-parameter  Margules  eqviations,  the  van  Laar . 
eqiiation  and  the  NRTL  (nonrandom,  two-liquid)  eqxaation  using  the 
IBM  360  computer.    The  results  obtained  were  compared  graphically 
on  plots  of  (1)  vapor  composition  versus  liquid  composition, 
(2)  liquid  phase  activity  coefficient  versus  liquid  composition 
and  (3)  total  vapor  pressure  versus  liquid  composition. 


